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Abstract 

CP-violating effects in the time-dependent angular distribution of the 
Bg — ► J/ip[—> £ + £~ ](/)[—> K + K~] decay products play a key role for the 
search of new physics. The hadronic Standard-Model uncertainties are 
related to doubly Cabibbo-suppressed penguin contributions and are 
usually assumed to be negligibly small. In view of recent results from 
the Tevatron and the quickly approaching start of the data taking 
at the LHC, we have a critical look at the impact of these terms, 
which could be enhanced through long-distance QCD phenomena, and 
explore the associated uncertainty for the measurement of the CP- 
violating B®-B® mixing phase. We point out that these effects can 
actually be controlled by means of an analysis of the time-dependent 
angular distribution of the B® — > J/ip[—> £ + £~]K*°[-^ ir + K~) decay 
products, and illustrate this through numerical studies. Moreover, we 
discuss SU (3)-breaking effects, which limit the theoretical accuracy of 
our method, and suggest internal consistency checks of SU(3). 
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1 Introduction 



The exploration of CP-violating effects in B s -meson decays offers a particu- 
larly promising probe for the search of New Physics (NP). In this respect, a 
key channel is B° s — > J/ip<f), which is the counterpart of the "golden" decay 
B% — > J I ipKs to measure the angle (3 in the unitarity triangle (UT) of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix. Since the B° s — > J/ip4> decay 
involves two vector mesons in the final state, the time-dependent angular 
distribution of the decay products of the vector mesons, J/tj) — > £ + £~ and 
<p — > K + K~, has to be measured in order to disentangle the admixture of 
different CP eigenstates [Hl2]. 

Within the Standard Model (SM), the CP- violating effects in the time- 
dependent B® — > J/'<P4 > angular distribution are expected to be small. On the 
other hand, a preferred mechanism to accommodate a measurement of non- 
vanishing CP asymmetries would be given by CP-violating NP contributions 
to B®-B® mixing (see, for instance, [3]). Recent results from the first tagged, 
time-dependent B° s — > J/ip<p analyses performed by the CDF [I] and D0 
[3] collaborations at the Tevatron (FNAL) may actually point towards this 
direction, and have led to quite some attention [B|. The B° s — > J/ip<fi decay 
is a main target of the LHCb experiment (CERN), which will soon start 
taking data and will allow us to explore the CP-violating phenomena in 
this transition with impressive accuracy [TJ: already with 2fb _1 of data, 
corresponding to one nominal year of operation, the experimental uncertainty 
for the B®-B® mixing phase S is expected to be a((p s ) cxp ~ 1°, and an 
upgrade of LHCb with an integrated luminosity of 100 fb _1 would eventually 
allow us to even reach a sensitivity of cr(0 s ) ex p ~ 0.2° [8]. 

In view of these exciting prospects, we have a closer look at the CP- 
violating effects in the time-dependent B® — > J/ip[—^ £ + £~](p{—>- K + K~] an- 
gular distribution that arise within the SM and limit the theoretical accuracy 
of the benchmark for the search for NP. Here the key role is played by pen- 
guin topologies, which are doubly Cabibbo suppressed and hence usually 
assumed to be negligible. However, these contributions cannot be calculated 
reliably from QCD, and could mimic CP-violating effects which might be 
misinterpreted as signals of NP in B®-B® mixing with a small but sizeable 
CP-violating NP phase. 

In the present paper, we point out that the penguin effects can actu- 
ally be controlled by means of an analysis of the angular distribution of 
B° s -> //^h £ + £-]K*°[^ n + K~] and its CP conjugate. Applying SU(3) 
flavour-symmetry arguments and neglecting penguin annihilation and ex- 
change topologies (which can be probed through B® — > J/ip<j)), the relevant 
hadronic parameters entering the B® — > J/ip<p observables can be deter- 
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mined, thereby allowing us to take them into account in the extraction of 
4> s - We suggest to perform a simultaneous analysis of the B° s — > J/ip4> and 
B® — > J/if)K*° channels at LHCb, and encourage the CDF and D0 collabo- 
rations to search for signals of this transition, as these would allow us to give 
first constraints on the penguin effects in B° s — > J/ipcf) and their impact on 
the extraction of the CP-violating B®-B° s mixing phase. Further information 
can be obtained from the B% — > J/ipp° decay, in particular for the resolution 
of a discrete ambiguity through experimental data. 

As pointed out in Ref. [5J, the data for CP violation in B® — > J/ijj7r 
and the branching ratio of this channel signal that such effects are size- 
able and soften the tension in the fit of the UT between its angle (3 and 
side Rb as determined through CP violation in B® — > J/ipKs^ decays and 
semileptonic b — > u,c transitions, respectively. In particular, the measure- 
ment of (3 has already reached a level of precision where subleading effects, 
i.e. doubly Cabibbo-suppressed penguin contributions, have to be included in 
order to match the experimental accuracy (see also Ref. [ID]). This feature 
strengthens the need to deal with such effects in analyses of CP violation 
in Bg — > J/ip<p as well. In particular, we expect that the penguin effects 
interfere constructively with mixing-induced CP violation and could lead to 
CP asymmetries as large as 0{— 10%), which would be significantly larger 
than the naive SM estimate of sin « —3% and could be well detected at 
LHCb. 

The outline of this paper is as follows: in Section [2j we give an overview of 
the B° s — > J/ip[— > £ + £"](f)[-^ K + K~] analysis and explore the impact of the 
penguin effects on the measurement of S , while we discuss the strategy to 
include the hadronic penguin contributions with the help of B° s — > J/i/jK*° 
in Section [3j This strategy is illustrated in Section [4} A detailed discussion 
of S"?7(3)-breaking effects and internal consistency checks that are offered by 
the observables of our decays into two vector mesons are given in Section [5} 
Finally, we summarize our conclusions in Section [6} 

2 Review of B° a J/V></> 

2.1 Structure of the Angular Distribution 

In contrast to the decay B® —>■ J/ipK§, we have to deal with two vector 
mesons in the final state of B® — > J/ip<f), which is an admixture of CP- 
odd and CP-even eigenstates. Using the angular distribution of the decay 
products of the vector mesons, the CP eigenstates can be disentangled. To 
this end, we introduce linear polarization states of the vector mesons, which 
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are longitudinal (0) or transverse to their directions of motion. In the latter 
case, the polarization states may be parallel (||) or perpendicular (X) to one 
another [llj. The time-dependent angular distribution of B® — ► J/ip<p takes 
the following general form [TJ: 

/(6,$,*;0 = X)» (fc) (©^.*)6 ( * ) (0, (!) 

k 

where the decay kinematics is described by the g^(Q, $, and the time- 
dependent coefficients b^ k \t) are given as 

^(t)! 2 (/e {0,11.x}), 

(2) 

Re{^(t)A||(t)}, lm{A%t)A ± (t)} (f G {0, ||}), 

with linear polarization amplitudes Ay = (( J / ijxt)) f\Hfjg\B® (t)), where TC e s 
is the relevant low-energy effective Hamiltonian. Here A±(t) describes a CP- 
odd final-state configuration, whereas A (t) and A»(t) correspond to CP-even 
final-state configurations. 

In the case of the CP-conjugate decay B° s — ► J/ip<f), we may write the 
angular distribution as 

7(9, $, t) = J2 & k \t)g( k \Q, $, (3) 

Since the meson content of the J/ip<f) state is the same whether it results 
from the 5° or B® decays, we may use the same angles 0, $ and \I> as in ([!]) 
to describe the kinematics of the decay products. Following these lines, the 
effects of CP transformations relating B° s — > (J/ip(j))f to B° s — > (J/ip<f))f are 
then taken into through the CP eigenvalues of the final-state configuration 
(J/ip(f>)f. Therefore the same functions g( fc )(0, $, are present in ([TJ and 
For the explicit form the of these quantities, see Ref. [TJ. 

2.2 Structure of the Decay Amplitudes 

As can be seen in Fig. [TJ colour-suppressed tree-diagram-like and penguin 
topologies contribute to the B° s — > J/ip<p decay within the SM. For a given 
final-state configuration / G {0, ||, X}, the B® — ► J/ip(f) decay amplitude can 
therefore be written as 

A(B° -> {J/W)f) = ^ 



Ac)f , 4(c)/ 

XI rp -|- XI p 



(*)/ 
p ) 



(4) 
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colour singlet 
exchange 




Figure 1: Decay topologies contributing to — > J /ipcf) in the SM. 



where the Xf = V js V* b are CKM factors, while A^ )f and A^ )f are CP- 
conserving strong amplitudes related to tree-diagram-like and penguin topolo- 
gies (with internal j G {u,c,t} quarks), respectively. Using the appropriate 
low-energy effective Hamiltonian, the latter quantities can be expressed in 
terms of linear combinations of perturbatively calculable Wilson coefficient 
functions and non-perturbative hadronic matrix elements of the correspond- 
ing four-quark operators, which are associated with large uncertainties. Using 
the CKM unitarity relation A^ 
we obtain 



— — Xu to eliminate the A^ factor, 



A(B° 8 -> (J/Vtf),) 



1-y )A, [l + eate»'e«] 



where 
and 



\ 2 A 



A ( " )f + A 



(c)f 



a f e ie f = R b 



A 



(«)/ 



A 



A 



(*)/ 



(*)/ 



A {c)f , Ac)f A (t)f 



(5) 



(6) 



(7) 



are CP-conserving hadronic parameters, while 

A = \V US \ = 0.22521 ± 0.00083, 
A=\V cb \/X 2 = 0.809 ±0.026, 



R b = (l-X 2 /2)\V ub /(XV cb )\ 
e = A 2 /(l - A 2 ) 



0.423±g;^| ± 0.029, 
0.053 



(9) 
(10) 
(11) 



are CKM parameters P, H2] , and the UT angle 7 flips its sign when consid- 
ering CP-conjugate processes: 

A(B° S -> ( J/W) t ) = rj,(l- y) A f I 1 + ^ df e- 1 ''} . 



(12) 
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Here 77/ is the CP eigenvalue of the final-state configuration (J/-00)/. 



2.3 Time-dependent Observables 



If we neglect CP violation in B° s -B° s oscillations, which can be probed through 
wrong-charge lepton asymmetries and is a tiny effect in the SM, the formalism 
of B®-B® mixing yields the following expressions [13J: 



r[f,t] = \A f (t)\ 2 + \A f (t)f 



_p(s), 
H 



\A f (t)\ 2 -\A f (t)\ 2 = 2e 



-r a t 



A f D cos(AM s t) + A f M sin(AM s t) 



(13) 
(14) 



where 1^ an d are the decay widths of the "light" and "heavy" B s mass 
eigenstates, respectively, Y s is their average, and AM S = - the 
difference of the mass eigenvalues. The labels "D" and "M" remind us that 
non- vanishing values of A^, and A^ are generated through direct and mixing- 
induced CP-violating effects, respectively 

Since the hadronic parameters cife t6f , which are essentially unknown, en- 
ter ^ and (12) in combination with the doubly Cabibbo-suppressed param- 
eter e, they are usually neglected. In this limit, we obtain 



T[f,t] = |A/}| 2 (l+r ?/ cos0 s )e- r L )i + (l-r/ / cos0 s )e- r H 



2 e rs *sim 



\A f (t)\ 2 - \A f (t)\ 2 = 2n f m 
where we have introduced the abbreviation J\ff 



,sin(AM s t), 



(15) 
(16) 



X 2 /2)A' f , and <j) s is 



the CP-violating B° S -B Q S mixing phase. In the ratio of the CP-violating rate 



difference (16), which requires the "tagging" of whether we had an initially, 



i.e. at time t — 0, present B° s or B® meson, and the "untagged" rate (15), 
the overall normalization |jV/| cancels, so that <p s can be extracted. For the 
corresponding time-dependences of the other observables provided by the 
angular distribution, see Ref. [2]. 

In the SM, we have </>f M = —2X 2 r] = —(2.12 ±0.11)°, where the numerical 
value follows from the current CKM fits [Hj. However, since B®-B® mixing 
is a strongly suppressed flavour-changing neutral-current (FCNC) process in 
the SM, it is a sensitive probe for NP effects in the TeV regime. Should 
new particles actually contribute to this phenomenon, the off-diagonal mass 
element of the mixing matrix is modified as follows [5]: 



M- 



12 



M 



s,SM 
12 



(1 + K.e"-) , 



(17) 
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Figure 2: The situation in the <t s -k s plane of the NP parameters for B° s -B° s 
mixing arising from current data and theory input, as discussed in the text. 

where k s measures the strength of the NP contribution with respect to the 
SM, and a s is a CP-violating NP phase. Consequently, we have 

AM S = AMf 1 1 + K s e ir7s | , (18) 

S = 0fl + NP = _ 2X 2 V + arg(1 + (lg) 

As discussed in Ref. [3], the values of p s = AM S / AMf M and can be 
converted into contours in the <j s -k s plane, which sets the parameter space 
for NP contributions to B®-B® mixing. 

For many years, only lower bounds on AM S were available from the LEP 
(CERN) experiments and SLD (SLAC) p3]. In 2006, the value of AM, could 
eventually be pinned down at the Tevatron [TB]. The current status can be 
summarized as follows: 

AM-! ( 1 8-56±0.87)ps- 1 (D0 collaboration [T7j), 

3 ~ \ (17.77 ± 0.10 ± O.O^ps- 1 (CDF collaboration [IE]). ^ } 

In order to determine the parameter p s from these measurements, the SM 
value of AM S is required, involving a hadronic parameter f^ s Bs s , which 
can be determined by means of lattice QCD techniques and introduces the 
corresponding uncertainties into the analysis. The HPQCD collaboration 
finds AMf M = 20.3(3.0) (0.8) ps" 1 P], which yields p s = 0.88 ± 0.13. 

Recently, following Refs. P Ej, the CDF [3] and D0 [5J collaborations 
have reported the first results from tagged, time-dependent analyses of the 
full three-angle distribution of the B° s — > J/if)[— > £ + £~]<p{— > K + K~] decay 
products. In an analysis by the UTfit collaboration [Bj, taking also other 
constraints into account, it is argued that these results may indicate CP- 
violating NP contributions to B®-B® mixing, which would immediately rule 
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out models with minimal flavour violation (MFV). Recently, a first average 
of the CDF and D0 data was presented by the Heavy Flavour Averaging 
Group (HFAG) [20J, corresponding to the following twofold solution: 

*„=(-44i£)° V (-1351S) . (21) 

In Fig. [2j we show - as an update of the analysis performed in Ref. [3] - the 
corresponding situation in the <j s -k s plane: the central hill-like region corre- 
sponds to p s , i.e. the mass difference AM S , while the two branches represent 
the twofold solution for S ; the overlap of the AM S and <p s constraints results 
in the two shaded allowed regions. It will be very interesting to monitor these 
measurements in the future. Fortunately, the B° s — > J analyses are very 
accessible at the LHCb experiment [7j, which will soon start taking data. 



2.4 Impact of Penguin Contributions 



The experimental results discussed in the previous section were obtained 
by assuming that the doubly Cabibbo-suppressed parameters ctfe l9f , which 
describe - sloppily speaking - the ratio of penguin to tree contributions, play 
a negligible role. In view of the search for NP signals, which requires a solid 
control of the SM effects, and the tremendous accuracy that can be achieved 
at LHCb, we generalize here the formulae to take also these contributions 
into account. 

Let us first have a look at the untagged observables. Following Ref. [13], 
we have 



R 



f 



|A/}| 2 (1+77/COS0, 



(22) 



+2ea/ cos Of {cos 7 + 77/ cos(0 s + 7)} + e a* {1 + 77/ cos(0 s + 27)} 



Rl 



so that 



= |A/>| 2 [(1 - 77/ COS0 S ) 

-2ea/ cos Of {cos 7 — 77/ cos(0 s + 7)} + e 2 aj {1 — r/f cos((p s + 27)} 



(23) 



r[/,t = 0] = Rl + R f u = 2|A/>| 2 1 + 2ea / cos^cos7 + e 2 aJ 

On the other hand, the CP-violating observables are given as follows: 

A£ = -2 \Nf\ 2 ea f sin 0/SU17, (25) 
= ^fWf? [sin0 s + 2ea / cos^ / sin(0 s + 7) + e 2 aJsin(0 s + 27)] .(26) 



(24) 
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Figure 3: Impact of the penguin parameter aj on the mixing-induced CP 
asymmetry r] f A f M for Of = 180° as function of the B°-B° s mixing phase (j) s . 



Note that (22) and (|23j) are not independent from (25) and (26), as 

, 2 



K) 2 + «)=^ H . 



(27) 



The ratio of the "tagged" rate difference (14) and the "untagged" rate 



(13) can be written as 



\A f (t)\ 2 - \ A f (t)\ 2 A f D cos(AM s t) + A f M sin(AM s t) 



\A f (t)\ 2 + \A f {t)\ 2 cosh(Ar s t/2) - ^ r sinh(Ar fl t/2) ' 



where AT S = T H 



W r( s ) 



, and 



R f u + R f L 



If we introduce 



Nf = 1 + 2ea/ cos 6^ cos7 + e a^- 



2 , r[/,t = o] 



2|W 

the corresponding observables take the following forms: 

— 2eaj-sin^sin7 



(28) 



(29) 



(30) 



(31) 
(32) 



Alt = [sin0 s + 2ea / cos0 / sin(0 s + 7) + e 2 a 2 f sm((j) s + 27)] , 

A^ Ar = — -rf- [cos S + 2eaf cosOf cos(<f) s + 7) + e 2 a 2 f cos(<p s + 27)] . (33) 
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Figure 4: Left panel: the dependence of A0{ on a/ for various values of Of] 
right panel: the dependence of on a/ for various values of Of. 

The measurement of A^ Ar relies on a sizeable value of the width difference 
AI\. Moreover, we have 



(AQ 2 + (A f M ) 2 + «) 2 = I- 



(34) 



For the extraction of S , the key observables are the A^; in Fig. |ij we 
illustrate the impact of the penguin parameter aj. Since a/e j6,/ is defined 
in ([7]) in such a way that Of is given by 180° if we assume factorization, we 
have used this value in order to calculate the curves shown in Fig. |3j For 
this strong phase the penguin effects are actually maximal in since only 
cos6f enters. On the other hand, the direct CP asymmetries would then 
vanish, as they are proportional to sinOf. 

We observe that in order to accommodate a value of (f) s ~ —44°, as given 
in (21), we would need a/ ~ 2.5-5, which appears completely unrealistic. 
However, since a,f suffers from large uncertainties, values as large as 0.5 ~ 1 
can a priori not be excluded. Should <p s take a value on the small side, 
these hadronic SM contributions would lead to a significant uncertainty in 
the extraction of the B®-B® mixing phase. 



In order to explore this effect in more detail, we use (31) and (32) to 
derive the following expression: 



1 - {M>) 



sin 



(35) 



where 



sm 



2ea/ cos 9 f sin 7 + e a j sin 27 



(36) 
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and 



so that 



l + 2ea /C os%cos7 + e 2 aJcos27 
cosA0^ = , (37) 



. 2eaf cos9f sin 7 + e 2 a 2 f sin 27 

tanA0{ = - f - f- 2 (38) 

1 + lea f cos v f cos 7 + e^aj cos 27 

It should be stressed that the shift A0{ of the B° s -B® mixing phase does not 
depend on the value of (j) s itself. In Fig. |4j we show the dependence of A</>{ 
on the penguin parameter aj for various values of 9f, and give - in order to 
monitor the corresponding direct CP asymmetries - a similar plot for A^. 
We observe that A0{ is of the same size as for af ~ 0.4, and that a 
value of df ~ 1 would induce a shift of A0{ ~ —5°. As can be seen in the 
left panel of Fig. 4 we have —0.05 <; ^ +0.05 for a/ ^ 1 and values of 
\9f — 180° I as large as 40°. Interestingly, as we expect cos9f < 0, the shift 
of <p s is expected to be negative as well, i.e. it would interfere constructively 
with 0^ M . These features are fully supported by our recent analysis of the 
B° — > J/ifjn° channel [9]. Consequently, it is important to get a handle on 
the penguin effects in the B° s — > J/ipcfi decay. 



3 The Control Channel B° a J/ipK* 

3.1 Structure of the Decay Amplitudes 

In Fig. [5J we show the decay topologies contributing to the B® — > J/ipK*° 
channel. The key difference with respect to the B° s — > J/i/xf) decay shown 
in Fig. [I] is that B® — > J/ipK*° is caused by b — > dec quark-level processes, 
whereas B° s — > J/i/j<fr originates from b — > sec transitions. Consequently, the 
CKM factors are different in these channels. In analogy to ([5]), we may write 

A(B° S -> (J/^K*°) f ) = XA' f \l - a' f e ie 'fe^] , (39) 

where A't and afre te f are the counterparts of the B° s — > {J/ip(f))f parameters 
introduced in (J6J) and 0, respectively. In contrast to (j5]), the latter param- 
eter does not enter p9~] ) in a doubly Cabibbo-suppressed way. Consequently, 
the B® — > J/if)K*° channel offers a sensitive probe for this quantity. If we 
apply the SU(3) flavour symmetry of strong interactions, we obtain 

\A f \ = \A' f \, (40) 
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c 




as well as 

a f = a' f , 6 f = 9' f . (41) 

In addition to 577(3) flavour-symmetry arguments we have here also assumed 
that penguin annihilation (PA) and exchange (E) topologies, which con- 
tribute to B® — > (J/ip<f>)f but have no counterpart in B° s — ► J/i^K* , play a 
negligible role. Fortunately, the importance of these topologies can be probed 
with the help of the B% — > (J/ip<f))f channel, which has amplitudes propor- 
tional to (PA + E)f. The Belle collaboration has recently reported the new 
upper bound of BR(B° d -> J/^0) < 9.4 x 1(T 7 (90% C.L.) [2TJ, which does 
not show any anomalous enhancement. The theoretical uncertainties associ- 
ated with the application of the 577(3) flavour symmetry will be discussed 
separately in Section [5j 



3.2 Observables 

In contrast to the B° s — > J/ip[-^ £ + £~}<f)[-^ K + K~] decay, the final states of 
B® — > J/if)[—> £ + £~]K*°[—> Tt + K~} and its CP conjugate are flavour-specific, 
i.e. the charges of the pions and kaons depend on whether we had a B® or 
B® meson in the initial state. Consequently, the time-dependent angular dis- 
tributions do not show CP violation due to interference between mixing and 
decay, i.e. the A^ observables introduced in (14) have no counterparts, and 
do not depend on the B° s -B° s mixing phase. However, untagged observables, 
as well as direct CP-violating asymmetries provide actually sufficient infor- 
mation to determine a'f and 9'j . In Appendix |Al we give the expressions for 
the time-dependent angular distributions, which allow the determination of 
the relevant observables. 
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Let us first discuss the untagged case, and introduce 



H f = 



A' 



mt = o]' 

r[/,* = o] 



1 — 2a'f cos Q'f cos 7 + a'j 
1 + 2ea/ cos Qf cos 7 + e 2 a 2 f ' 



(42) 



where T[f,t = Of is the B° s -> J/if)K*° counterpart of Q. Using Q, we 
may extract if/ from the untagged observables. Moreover, using also (41 ), we 
can determine a'f as a function of Q'f with the help of the following formulae: 



u H <±Ju 2 H -v, 



H 



/' 



where 



and 



U 



1 + eH 



l H 



cos 0'f cos 7, 



V Hf 



l-H 



f 



(43) 



(44) 



(45) 



1 - e>H f 

Here the main uncertainty is associated with the determination of Hf, which 



relies on (40). In Subsection 5.2, we have a closer look at the corresponding 



Si7(3)-breaking corrections, and give numerical results for the extraction of 
the Hf from the untagged observables. On the other hand, thanks to the e 



terms in (42), the impact of corrections to (41) is tiny. 



Another useful quantity is offered by the direct CP asymmetry 

2a' f sin Q' f sin 7 



If 



2a! f cos Q'f cos 7 + a'i 



12 > 



(46) 



which can be extracted from a rate difference; it takes the same form as (28) 
for t = 0. In analogy to Hf, also allows us to determine function 
of Q'f. To this end, we may again use (43), with the following replacements: 



U. 



H, 



u Af > 



cos Q'f cos 7 + 



sin Q'f sin 7 



AC 



v, 



V, 



1. 



(47) 



It should be emphasized that the corresponding curve in the Q'f-a!f plane is 



theoretically clean, whereas that described by (43) is affected in particular 



by the S , f/(3)-breaking effects entering the determination of Hf. 

The intersection of the Hf and A(y contours allows us then to extract 



and Q'r from the data. Finally, applying (41) and the results derived in 



Section 2.4 , we can include the penguin effects in the determination of the 
B®-B® mixing phase. Let us first illustrate this method in the next section 
by discussing a numerical example before giving a detailed discussion of the 
relevant SU (3)-breaking effects in Section [5] 
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Figure 6: Illustration of the contours in the d'f-off plane for an example, as 
discussed in the text. 



4 A Numerical Example 

For the illustration of the strategy discussed above, we assume 7 = 65°, 
and hadronic parameters given by a'* = 0.4 and 9'^ = 220°, yielding the 

observables Hf = 1.44 and = —0.33. These input values are consistent 
with the ranges of the 5° — > J/iprr parameters a' £ [0.15,0.67] and 6' £ 
[174°, 213°] found in Ref. j9]; we expect a picture for a'j and 9'j that is similar 
to the one for their B® — > J/tpir counterparts. 

In Fig. [6j we show the contours in the O'j-a'j plane that arise in this ex- 
ample. We observe that a twofold solution emerges for (6'f,a'f), which can 
be resolved through the sign of cos#j. Theoretically, we expect a negative 
value of this quantity, which is also supported by the B® — > J/ipn° data. In 
order to resolve this ambiguity experimentally, we need an additional observ- 
able, which would be provided by mixing-induced CP violation. Since the 
B® — > J/ip[-^ £ + £~]K*°[-^ n + K~] processes have flavour-specific final states, 
they do not show this phenomenon. On the other hand, mixing-induced CP 
violation would arise in B° s — > J/ip{—>- £ + £~]K*°[— > 7r°Ks^} modes, in anal- 
ogy to B° — > J/ip[-^ £ + £~]K*°[— > 7i°K SjL ] processes [I]. Unfortunately, it is 
essentially impossible to study the corresponding experimental signatures in 
a hadronic environment, i.e. at the Tevatron or LHC. 

However, we may alternatively use the B^ — > J/ipp° channel [13] . which 
can be obtained from B® — > J/ipK*° by replacing the strange spectator quark 
through a down quark, as can be seen in Fig. [5j In this case, the final state 
is an admixture of different CP eigenstates, in analogy to B® — > J/ip(fi, and 
we can extract the following mixing-induced CP asymmetry from the time- 
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Figure 7: Illustration of the resolution of the twofold ambiguity in Fig. [6] 
through the mixing-induced CP violation in B® — ► J/ipp°. 



dependent angular distribution: 
+77/ 



sin <2> d 



2 a 'f cos 6*j sin( 



+ 7) + a'* sin(0 d 



+ 2 7 ) 



1 — 2a! * cos #y cos 7 + a'j 



(48) 



where 77/ is the CP eigenvalue of the final-state configuration /, i.e., 7/0, 
m = +1 and r/± = —1, whereas <pd = (42.4j^ 7 )° denotes the B^-B® mixing 
phase [9]; for simplicity, we have also denoted the B® — > J/ipp° hadronic 
parameters by a'* and 9'f, as we expect them to be approximately equal to 
those of B® — > J/ipK*° thanks to the SU(3) flavour symmetry. Using (43) 
with the replacements 



U, 



u Af > = 



sin 



M COS 7 



sm{(p d 



+ 2 7 ) 



M 



sin <p d 



- A f ' 

Af 



cos 9' , 



(49) 



(50) 



sin(0 d + 2 7 ) -A£ 

the measurement of the mixing-induced CP asymmetry A^ allows us to fix 
another contour in the 9'f-a'f plane. If we consider the example given above 

with <pd = 42.4°, we obtain r/jA^ = 0.90, which results in the contours 
shown in Fig. [7j We see that the twofold ambiguity in the determination of 
the hadronic parameters can now be resolved, thereby leaving us with our 
input values. 

Since the width of the p° is three-times larger than that of the K*°, 
the B° s — > J/ipK*° control channel should be experimentally better ac- 
cessible than B® — > J/ipp°. Moreover, if we neglect S , t/(3)-breaking ef- 
fects due to the different spectator quarks, we expect the simple relation 
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Figure 8: Situation in the A0{ plane for the example in Figs. [6]and[7j 



BR(£° -> J/ijK*°) ~2x BR(B° d -> J/^p°) = (4.6 ± 0.4) x 10" 5 [20]. How- 
ever, already a rather crude measurement of the mixing-induced CP- violating 
J/ipp would be sufficient to resolve the ambiguity in the 
In particular, the expected negative value of cos 9'j- 
that are larger than sin0 rf = 0.67. 
Such a pattern emerges actually in the measurement of the mixing-induced 



observables of B d 
extraction of a'* and 9'j 

would be indicated by values of f]fA J M that are larger than sin 



CP violation of B° d 



In Fig. [8 

by means o 



we convert the contours in the O'f-a'f plane into the #/-A0{ space 



(36)-(38) and (41). We observe that, in this specific example, 
the shift of the B® -B® mixing phase through the penguin effects is given by 
A0{ = — 1.7°. If we assume the SM, the mixing-induced CP asymmetries of 
B° — > J ft/Hp represented by (35) would be given by r/fA^ = —6.7%, which 
is about twice as large as the SM value. At LHCb, such CP asymmetries 
could be detected with about 4 a significance after collecting 2 fb _1 of data, 
corresponding to one nominal year of operation, and with about 20 a at an 
upgrade of this experiment with 100fb _1 integrated luminosity. However, 
without the control of the hadronic penguin effects through a simultaneous 
analysis of the B° s — > J/ipK* channel as proposed above, this result would 
be misinterpreted cLS cL SI gnal of physics beyond the SM. In this context it 
is important to emphasize that we expect and A0{ to have the same 
negative sign, thereby leading to constructive interference. In the opposite 
case, i.e. with a positive value of A(f>{, the SM picture of expecting vanish- 
ingly small CP violation in B® — > J/ip<p would be much more robust with 
respect to the hadronic penguin uncertainties. It cannot be excluded that 
the hadronic penguin effects are actually more significant than in our exam- 
ple, and could lead to 77/J% ~ —10%. This feature is fully supported by the 
picture emerging from the current B® — > J/ipn data [9]. 
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Figure 9: The maximal shift of <p s as function of Hf for various values of 7. 

In view of these findings, it would be very desirable to search for the 
B® — > J/ipK*° decay at the Tevatron. Already information on Hf would 
allow us to put first valuable constraints on the shift A0{ . As we have shown 
in Fig. |9j these observables will put a first upper bound on A0{. Once direct 
CP violation in the B° s — > J/ipK*° angular distribution is measured, A<p{ 
can be fully pinned down, as we have shown above. 



5 >S't/(3)-Breaking Effects 

5.1 General Remarks and lj— <fi Mixing 

The main theoretical uncertainty of the strategy proposed in the present pa- 



per is related to S'?7(3)-breaking effects which affect the relations in (40) and 



(41). The following discussion of SU(3) breaking is parallel to our previous 
investigation [5]; however, here we deal with vector meson final states. On 
the one hand this simplifies the discussion since S'[/(3)-breaking effects in 
the vector meson octett seem to be smaller; on the other hand, we have a <fi 
state which is believed to be to a good approxmation an ss state and hence 
a superposition of SU(3) eigenstates. 

The SU(3) nonet contains three neutral, non-strange states. Assuming 
isospin to be a good symmetry, one of these states is the neutral I = 1 p 
meson with the quark decomposition po — (uu — dd) / y/2. The other two 
states are isosinglets and are given by 0o — («h + dd + ss)/y/3, which is an 
SU(3) singlet, and <ps = {uu + dd — 2ss)/4 belonging to the SU(3) octet. In 
case of unbroken SU (3) symmetry, the amplitudes for processes involving the 
members of the octet are related to one another, while the singlet remains 
separate. 
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However, there is good indiction that the physical state is to a good 
approximation a pure ss state and hence a superposition of the singlet O 
and the octet 8 . The orthogonal state cu = [uu + dd) /y/2 is the isoscalar 
lo meson. Up to small mixing between u and 0, which has been discussed 
recently in the context of B decays in Ref. [22], and which is too small to be 
relevant here, these are the (strong) mass eigenstates. 

The way we apply the SU (3) symmetry is to assume that the matrix 
elements of = ss are related to the corresponding matrix elements of the 
members of the octet. In other words, we shall assume that the form factors 
of the Bg — > transition are the same as the ones for the B® —>■ K*° decay. 
Strictly speaking, this goes beyond the SU(3) symmetry assumption since 
we relate octet and singlet components. Lacking any detailed information 
on the quality of such an asumption we have to rely, e.g., on QCD sum rule 
estimates which indicate that the strong dynamics in the = (ss)s=i state 
are very similar to K*° = (ds)s=i', in fact, we shall rely on QCD sum rules in 
Section ^2 to discuss the deviations from our assumption. In this context, it 
should be emphasized again that we have also to neglect penguin annihilation 
and exchange topologies, which can be probed through the B® — > J/ip4> decay. 

In a recent paper |23J, it is argued in detail that the relations between 
Bg — > J/ip4> and B® — > J/ipK*° following from flavour symmetry [24J are 
likely to be quite reliable, so that using strong phase information from the 
B® — > J/ipK*° channel in the analysis of B° s — > J/ip(j) ls justified. Also here 
we have to assume that the matrix elements of = ss are related to the 
corresponding matrix elements of the members of the octet. 

It is very difficult to get a reliable estimate of the SU(3) breaking for the 
non-leptonic decays at hand. It is known from the corresponding processes 
with pseudoscalar final states that the decays with J/ip in the final state 
are dominated by non-factorizable contributions; it is not even clear how to 
factorize the penguin contributions in the decays we are considering. How- 
ever, in the case of B^ — > tttt, ttK, KK decays, we encounter sizeable non- 
factorizable effects, whereas the data do not indicate large S'[/(3)-breaking 
effects of this kind |25j. In particular, considering the counterpart of the Hf 
quantities introduced in the present paper for the B® — > K + K~ , B® — > 7r + 7r~ 
system, a calculation of the relevant form-factor ratio by means of QCD sum 
rule techniques [26J yields good agreement with the current data that would 
be spoiled by large non-factorizable, S"?7(3)-breaking effects. 

This empirical behaviour gives us conficence that our estimate of the 
S"?7(3)-breaking effects for the extraction of the Hf from the data given in 
the next subsection, which relies on a QCD sum rule analysis of the relevant 
form factors as well, describes the leading corrections. 
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5.2 SU(3) Breaking in the Extraction of Hf 

In order to calculate the S77(3)-breaking corrections to the amplitude ra- 
tios \A.f/A!f \ that are required for the extraction of the Hf from the data 
(see (42)), we apply the formulae given in Ref. pp. The linear polarization 



amplitudes of the B® — > J/ip4> channel at time t = can be written as 

Aq(0) = -xa - (x 2 - l)b 

A||(0) = V^o (51) 
A ± {0) = ^/2(x 2 - 1) c 

with 

X - mj/ ^- 2m m mt ' (52) 
where the "factorized" contributions are given by 

&fact = ^|A« (Cf (ai) + CfM)B**, (53) 

Cfact = ^|A( s )(Cf( Ai )+Cf( / i))Cr. 

Here we have neglected the doubly Cabibbo-suppressed penguin corrections, 
as our target are the overall amplitudes Af] Gp is Fermi's constant, Ac the 
CKM factor introduced after Q, and the Cf ff (/i) are the "effective" Wilson 
coefficient functions introduced in Ref. pp. Moreover, we have 

Af ct = -fj,^mj^{m Bs + m^)A^{m 2 ^) , 



Bf* = 2 J Z XL * Af\m 2 m ) , (54) 



fj/^ 2 j/^<f> 
m Bs + 



where Af^g 2 ) and V Ss ^((7 2 ) are the form factors of the quark-current ma- 
trix elements of the B s — > <p transition, with q denoting the momentum 
transferred by the quark current. In the case of the B° s — > J/ipK*° channel, 
we need correspondingly the B s — > K*° transition form factors, and have to 



replace — > K*° in (52) and (54). 

An analysis of these form factors was performed in Ref. |27j. Light cone 
QCD sum rules allow an estimate of the values of the form factors at q 2 = 
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v = 4> 


V = K* 


A"'~* v (m 2 ji^) 
A B ^ V {m% f ) 


0.42 ± 0.06 
0.38 ±0.06 
0.82 ±0.12 


0.33 ± 0.05 
0.32 ±0.05 
0.62 ± 0.09 



Table 1 : Collection of the relevant B s — > V form factors at q 2 
the results of Ref. [27J and assuming an uncertainty of 15%. 



0. In order to obtain the value of the form factor at a different q 2 , such 

J2. 



as q 



m 



as in (54), we have to make an extrapolation using some 



parametrization of the rorm factor. If we use the functional forms suggested 

in Ref. [27] and assume an uncertainty of 15%, we obtain the form factors at 
2 ° 

r 



m 



collected in Table [TJ and the following SU (3)-breaking ratios: 





->K* 


( m2 J/i>) _ 




-*( 




Af"' 




_ 


A*' 


-+( 




V B 8 - 


->K* 






-+( 





0.78 ±0.08 



0.84 ±0.07 



0.76 ±0.15 . 



(55) 



Using then (51), we obtain the following numerical results, which allow the 



extraction of the Hf from the untagged rates with the help of (42) 

2 

= 0.42 ±0.27 , 



A' 



Ao 



A\\ 
A\ 



Ai 



0.70 ±0.29 



0.38 ±0.16 



(56) 



Note that in order to calculate |^4q/^4 | 2 ) we need the 
factors given in Table [TJ 



form 



5.3 SU(3) Breaking in a' f = a f and Q' f = f 

If we use the B® 



J/ipK*° observables as discussed in Section [3j we can 

J/ip4> counterparts af and 



extract a'j and 6'f from the data. Since their B Q S 
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Of enter in Hf in combination with the tiny parameter e, this determination 



is essentially unaffected by corrections to (41); the main corrections enter 
through the value of Hf, which requires the amplitude ratios [Af/A'A, with 
the £?7(3)-breaking corrections estimated in the previous subsection. 



When calculating the shifts A0{, we have to use the relations in (41). 
However, one has to keep in mind that sizable non-factorizable effects could 
induce £?7(3)-breaking corrections. Their impact on the determination of 



A0{ can be easily inferred from (38). Neglecting terms of order e 2 , we have 
a linear dependence on cifCosOf. Consequently, corrections to the left-hand 



side of (41) propagate linearly, while ££/(3)-breaking effects in the strong 
phases will generally lead to an asymmetric uncertainty for A<p{. 

In the analysis of the B® — > J/ip7r° data in Ref. [9], the impact of SU(3)- 
breaking corrections was explored by setting a = C,a' and uncorrelating the 
strong phases 9 and 9' of the — > J/ipK° and B® — > J/i/jtt decays, respec- 
tively. Even when allowing for £ G [0.5, 1.5] and 9,9' G [90,270]° in the cor- 
responding fit, and using a 50% increased error for the relevant form-factor 
ratio to explore the impact of dramatic non-factorizable, ££7(3)-breaking 
contributions to |^4/^'|, the picture emerging from the global fit is not sig- 
nificantly changed. To be specific, A0^ G [—6.7,0.0]° arises when allowing 
for such large ££7(3)-breaking corrections, whereas Aifid G [—3.9,-0.8]° in 
the case with £ = 1 and 9 = 9'. We expect a similar situation for A</>{. 

5.4 Internal Consistency Checks of SU (3) 

The advantage of B decays into two vector mesons is that many more ob- 
servables are offered by the angular distribution of their decay products than 
in the case of B — > PP or B — > PV decays (P and V denote generically 
pseudoscalar and vector mesons, respectively). This comment applies also to 
the decays considered in the present paper, and allows us to perform internal 
consistency checks of the SU(3) flavour symmetry. 

A very first internal test follows from a comparison of the different values 
of the B° s -B° s mixing phase (f) s following from the three polarization states 
/ G {0, ||, _L}. Obviously, these values should agree with one another. In fact, 
even more quantitative tests of £77(3) breaking can be performed. The point 
is that we may choose one of the three linear polarization states to extract 



4> s from (35), taking the shift A0{ through the penguin effects into account. 
Using then the B® J /ip4> observables and of the remaining two 
polarization states, the knowledge of <p s allows us to extract the corresponding 



shifts A0{ from (35). With the help of (38), we can then convert the values 



of the A0{ into contours in the Of— (if plane. To this end, we have simply to 
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Figure 10: The extraction of aj and Of from A<j){ and the direct CP asym- 
metry A(^ for internal consistency checks of the SU(3) flavour symmetry, as 
described in the text. The example corresponds to aj = 0.4 and Of = 220° 
with 7 = 65°, as in the previous numerical illustrations. 



make the following replacements in (43) 



U 



_ { sin 7 — cos 7 tan A0{ 
)S 27 tan A<p{ — sin 27 , 



cos 



V, 



V A<f>l 



tan A0{ 



(57) 



(5f 



as well as a'f 
substitution, 0', 



f 



cos 27 tan A<p s — sin 27 
ecif. Moreover, if we replace, in addition to the latter 
> 180° + f and A£ — > in (47), the direct CP asymmetry 



in Bg — > J/ip(j) can be converted into a contour in the Of-af plane as well 



It should be stressed that these constructions are valid exactly. In Fig. 10 



we illustrate how this works by considering again the numerical example 
specified in Section |4j 

The values of the hadronic B° s — > J/ip4> parameters af and Of allow us 
then to perform an internal consistency check of the SU(3) flavour symmetry 
by comparing with the values of a'f and O'f following from the B® — > J/ipK*° 
strategy proposed in Section [3] Another test is offered by the following 
relations: 

A f D = -eH f A f D , (59) 



which rely on (41). Needless to note, the practical usefulness of these consis- 



tency checks depends on the values of the observables that will eventually be 
measured by LHCb. We strongly encourage detailed feasibility studies and 
look forward to confronting these considerations with real data soon. 
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6 Conclusions 



Studies of CP-violating effects in the time-dependent angular distribution of 
B® — > J/^[— > > processes have recently received consider- 

able attention in view of first tagged measurements at the Tevatron, and are a 
central target of the LHCb experiment which will soon start taking data. We 
have pointed out that hadronic effects, which are due to doubly Cabibbo- 
suppressed penguin contributions that are usually neglected, could induce 
mixing-induced CP-violating effects as large as 0(— 10%). Without the con- 
trol of these penguin contributions, which cannot be calcuated reliably from 
QCD, such CP-violating effects, which can be detected with excellent signif- 
icance by LHCb, would be misinterpreted as CP- violating NP contributions 
to B° s -B° s mixing. 

In the present paper, we have proposed a strategy to include these con- 
tributions with the help of a measurement of the angular distribution of the 
B° s — > J/il)[—> £ + £~]K*°[-^ 7i + K~] decay products, and have illustrated this 
by means of a numerical example. We strongly suggest a search for this 
control channel at the Tevatron in order to obtain first constraints on the 
penguin effects in the B° s — > J/ip<fi analysis. The tremendous accuracy that 
can be achieved at LHCb and a possible future upgrade of this experiment 
makes it mandatory to include these penguin contributions. 

Acknowledgements: S.F. and T.M. acknowledge the support by the Ger- 
man Ministry of Research (BMBF, Contract No. 05HT6PSA), and R.F. 
would like to thank Guido Martinelli for the hospitality at the Universita 
di Roma "La Sapienza" . 
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Appendix 

A Time-dependent Angular Distributions of 
B° s J/i/jK*° and CP Conjugates 

Following Ref. [T], we introduce the following set of trigonometric functions: 

fx — 2 cos 2 t/> (1 — sin 2 cos 2 <£>) 

/2 = sin 2 t/i (1 — sin 2 6 1 sin 2 y?) 

/ 3 = sin 2 ip sin 2 5 

/4 = sin 2 ^ sin 29 sin </? ^ ' 

fs = (1/V2) sin 2^ sm 2 #sin2y2 

/ 6 = (l/y/2) sin 2^ sin 2$ cos p. 

If we use the notation Af = A(B° S — > (J/ipK*°) f) for the unevolved amplitude 



in (39) and Af for its CP conjugate, we obtain 



d 3 T[B°Jt) -> JM(-> vr+if-)] 9 r 1/A „ /n , /asj .m rt 

1 sKJ > K >± = — cosh AT a t/2) + cos AM s t e" 1 ^ 

a cos 6 dip d cos ip b4n 

x {Al^ol 2 + / 2 |A|| | 2 + f 3 \A ± \ 2 - him + / 5 Re + him (A* A ± )} 

(61) 

^ l a = cosh Ar s t/2 + cos(AM s t) e"^ 

dcosu dp d cos ip oAn 

x {Al^ol 2 + i^iil 2 + y*3 1 1 2 + him (A*\A±) + f 5 Re(AlA\\) - him (A* A ± )} 

(62) 

d 3 T[B°Jt) -> JM(-> £ + £~)K*°(^ 7T-K+)] 9 r 1/AT , /0 , , A „, X1 r * 

3 fl w 3 / = cosh Ar^/2 - cos AM s t e"^ 

a cos dp d cos ip 64n 

x {/ll^ol 2 + /2|-4||| 2 + h\A±\ 2 + him (AjjAi) + / 5 Re(AJi4||) - / 6 Im 

(63) 

, „ , — , , = [cosh(Ar s t/2) - cos(AM s t) e~ Vat 

a cost* dp d cos ip 647r 

x {h\A \ 2 + f 2 \A {l | 2 + f 3 \A ± \ 2 - h lm (AfAJ + / 5 Re + / 6 Im (A* A ± )} . 

(64) 
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In the case of Ar s — > 0, we have 

cosh(Ar s t/2) + cos(AM s t) -> 2 cos 2 (AM s t/2), (65) 

cosh(Ar s t/2) - cos(AM s t) -> 2 sin 2 (AM s t/2). (66) 

Consequently, the expressions listed above reduce to those given in Ref. pQ 
for the flavour-specific — > J/ip{—> £ + £~]K*[—+ K -k^\ modes with the 
assumption of \Af \ = \Af\. 
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